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Abstract 

Background: The study of synchronization among genetic oscillators is essential for the understanding 
of the rhythmic phenomena of living organisms at both molecular and cellular levels. Genetic networks are 
intrinsically noisy due to natural random intra- and inter-cellular fluctuations. Therefore, it is important to 
study the effects of noise perturbation on the synchronous dynamics of genetic oscillators. From the synthetic 
biology viewpoint, it is also important to implement biological systems that minimizing the negative influence 
of the perturbations. 

Results: In this paper, based on systems biology approach, we provide a general theoretical result on 
the synchronization of genetic oscillators with stochastic perturbations. By exploiting the specific properties 
of many genetic oscillator models, we provide an easy-verified sufficient condition for the stochastic syn- 
chronization of coupled genetic oscillators, based on the Lur'e system approach in control theory. A design 
principle for minimizing the influence of noise is also presented. To demonstrate the effectiveness of our 
theoretical results, a population of coupled repressillators is adopted as a numerical example. 

Conclusion: In summary, we present an efficient theoretical method for analyzing the synchronization 
of genetic oscillator networks, which is helpful for understanding and testing the synchronization phenomena 
in biological organisms. Besides, the results are actually applicable to general oscillator networks. 
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Background 

Elucidating the collective dynamics of coupled genetic oscillators not only is important for the understanding of 
the rhythmic phenomena of living organisms, but also has many potential applications in bioengineering areas. 
So far, many researchers have studied the synchronization in genetic networks from the aspects of experiment, 
numerical simulation and theoretical analysis. For instance, in [T], the authors experimentally investigated the 
synchronization of cellular clock in the suprachiasmatic nucleus (SCN); in [USE], the synchronization are 
studied in biological networks of identical genetic oscillators; and in [5j O [7], the synchronization for coupled 
nonidcntical genetic oscillators is investigated. Gene regulation is an intrinsically noisy process, which is subject 
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to intracellular and extracellular noise perturbations and environment fluctuations [H O EH H3 HH [Hj ■ Such 
cellular noises will undoubtedly affect the dynamics of the networks both quantitatively and qualitatively. In 13J, 
the authors numerically studied the cooperative behaviors of a multicell system with noise perturbations. But 
to our knowledge, the synchronization properties of stochastic genetic networks have not yet been theoretically 
studied. 

This paper aims to provide a theoretical result for the synchronization of coupled genetic oscillators with 
noise perturbations, based on control theory approach. We first provide a general theoretical result for the 
stochastic synchronization of coupled oscillators. After that, by taking the specific structure of many model 
genetic oscillators into account, we present a sufficient condition for the stochastic synchronization in terms of 
linear matrix inequalities (LMIs) |f 5j . which are very easy to be verified numerically. To our knowledge, the 
synchronization of complex oscillator networks with noise perturbations, even not in the biological context, has 
not yet been fully studied. Recently, it was found that many biological networks are complex networks with 
small- world and scale-free properties p~6l [UJ . Our method is also applicable to genetic oscillator networks with 
complex topology, directed and weighted couplings. To demonstrate the effectiveness of the theoretical results, 
we present a simulation example of coupled repressilators. Throughout this paper, matrix U € j{ NxN j s defined 
as an irreducible matrices with zero row sums, whose off-diagonal elements are all non-positive, and the other 
notations are defined in the Appendix A. 



Results 

Theoretical Results 

Since we know very little about how the cellular noises act on the genetic networks, a simple way to incorporate 
random effects is to assume that certain noises randomly perturb the genetic networks in an additive manner. 
We consider the following networks of N coupled genetic oscillators with random noise perturbations 

dx (t) N 

-jp = F{ Xi (t)) + y £ / G ij Dx j {t) + Vi(t)m(t), i =!,■■■ ,N 

3=1 

where F(.) defines the dynamics of individual oscillators, Vi(t) 6 R nxl is called the noise intensity vector, 
belongs to L2[0, oo). As we will see in the following analysis, the results hold no matter what Vi(t) is and 
no matter where it is introduced, so we do not explicitly express the form of Vi(t) here. Vi(t) can also be a 
function of the variables (if so, some minor modifications are needed in the following). rii(t) is a scalar zero 
mean Gaussian white noise process. Recall that the time derivative of a Wiener process is a white noise process 
19J, hence we can define dwi(t) = rii(t)dt, where Wi(t) is a scalar Wiener process. Thus, the above equation 
can be rewritten as the following stochastic differential equation form: 



dxi(t) 



x 

3 = 1 



dt + Vi(t)dwi(t), i = 1, • • • , N. (1) 



The work can be easily extended to the case that Vi(t) £ Ji nxl i arLC j m(t) — [nn(t),- ■ ■ ) nu i (t)] T be an U- 
dimensional mutually independent zero mean Gaussian white noise process. D G Ji nxn defines the coupling 
between two genetic oscillators. G = (Gij)NxN is the coupling matrix of the network. If there is a link from 
oscillator j to oscillator i (j ^ i), then Gij equals to a positive constant denoting the coupling strength of this 
link; otherwise, Gij — 0; Gu — — Y^jLx j=a Gij ■ Matrix G defines the coupling topology, direction, and the 
coupling strength of the network. 
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For network (1), a natural attempt is to study the mean-square asymptotic synchronization. But existing 
experimental results show that usually the genetic oscillators can not achieve mean-square synchronization 
(see, e.g. experimental results in [1] and Appendix B for a theoretical discussion). Analogue to the stochastic 
stability with disturbance attenuation (see, e.g. [20] ) , we give a less restrictive (but more realistic) definition of 
the stochastic synchronization as follows: 

Definition 1: For a given scalar 7 > 0, the network (1) is said to be stochastically synchronous (under 
the combination matrix U) with disturbance attenuation 7 if the network without disturbance (vi — 0, Vi) is 
asymptotically synchronous, and under the same initial conditions for all oscillators, 

£(-[/^IN(t)-^(i)||| 2 <7Z>WH2' W 

i<j i 

for all nonzero Vi(t), where || • |||. = E (J °° | • \ 2 dt). Here, by introducing the combination matrix U, we can 
flexibly select the form of the matrix to obtain different error combinations. 

By using the techniques described in the Appendix C, we know that if there exist matrices P > 0, T E R nxn 
and U, and a scalar p > 0, such that the following conditions are satisfied, 

S 2 = 2( yi - y 2 ) T P[F( yi ) - F(y 2 ) - T{y x - y 2 )} + L{y x - y 2 ) T ( yi - y 2 ) < 0, V Vl ,y 2 G R n ( Vl ± y 2 ), 

(U ® P)(G ®D + I<g>T) + (G®D + I® T) T (U ® P) < 0, (3) 
U®P <pl, 

then, the network (1) will achieve stochastic synchronization with disturbance attenuation 7. 

The above condition (3) is a general result for the stochastic synchronization of coupled oscillators. But we do 
not have general efficient method for verifying the first inequality in (3) for arbitrary F due to its nonlinearity. 
Next we consider a special structure of genetic oscillators to obtain an easy- verified result. 

Genetic oscillators are biochemically dynamical networks, which can usually be modelled as nonlinear dynam- 
ical systems. Mathematically many genetic oscillators can be expressed in the form of multiple additive terms, 
and the terms are monotonic functions of each variable, which particularly, are of linear, Michaelis-Menten and 
Hill forms. In our previous papers [7, 23 a , we have taken such special structure properties of gene networks into 
account, and have shown that these genetic oscillators can be transformed into Lur'e form and can be further 
analyzed by using Lur'e system method in control theory |22j . In this paper, we also consider such special 
structure. To make our paper more understandable and self-contained, we will first introduce the approach 
briefly, and after that we will analyze the stochastic genetic oscillator networks theoretically. We consider the 
following general form of genetic oscillator: 

ij{t) = Ay{t) + B l f 1 {y{t)) + B 2 f 2 {y{t)), (4) 

where y(t) £ R n represents the concentrations of proteins, RNAs and chemical complexes, A, B>i and B 2 are 
matrices in R nxn , fi(y(t)) = [fn(yi(t)), ■ ■ ■ , fin(yn(t))] T with fij(yj(t)) as a monotonic increasing function 
of the form fij( yj (t)) = { Vj (t) / /[I + (y 3 (t)/ (3 xj ) H ^], and f 2 (y(t)) = [f 21 ( Vl (i)), ■ • ■ , f 2n (y n (t))] T with 

f2j{yj{t)) as a monotonic decreasing function of the form f 2 j(yj(tj) = 1/[1 + (yj(t)/ f3 2 j) H2j ], where H±j and 
H 2 j are the Hill coefficients. Genetic oscillators of the form (4) is by no mean peculiar. Many well-known 
genetic system models can be represented in this form, such as the Goodwin model [24], the repressilator [25] . 
the toggle switch [26] , and the circadian oscillators [27] . Undoubtedly, this work can be easily generalized to the 
case of y(t) = Ay(t) + X)j=i Bjfj(y(t)), where there are more than two nonlinear terms in each equation of the 
genetic oscillator. From the synthetic biology viewpoint, genetic oscillators with only linear, Michaelis-Menten 
and Hill terms can also be implemented experimentally. 
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To avoid confusion, we let the jth column of Bi 2 be zeros if fij, 2 j = 0. Since f2j{Vj{t)) 



1 - = 1 ~ 9j(Vi(t)\ and letting /(.) = /i(.), we can rewrite (4) as follows: 

y(t) = Ay(t) + B 1 f(y(t)) - B 2 g{y{t)) + B 2 . (5) 
Obviously, fi and gi satisfy the sector conditions: < ^'( a ^~^'( b ) < ku, < 9i ^^{b) <- ^ 2i) Qr e q U i va i en tly, 

(/<(<*) - fi(b)Wi(a) - Mb)) - k u (a - b)] < 0, 

( 9i (a) - &(&))[(& (a) - 9l (b)) - k 2l {a - b)} < 0, (6) 
Va, b € R (a ^ b); % = 1, • • • , n, 

Recall that a Lur'e system is a linear dynamic system, feedback interconnected to a static nonlinearity that 
satisfies a sector condition [22]. Hence, the genetic oscillator (5) can be seen as a Lur'e system, which can be 
investigated by using the fruitful Lur'e system approach in control theory. 

By substituting the individual genetic oscillator dynamics (5) for F in the network (1), we obtain the following 
network of ./V coupled genetic oscillators: 

dxi(t) = [Axi(t) + Bif(xi(t)) - B 2 g{x l {t)) + B 2 + Eli G lJ Dx 3 {t)]dt + Vi(t)dwi(t), 

(7) 

i = !,-■■ ,N. 
For this network, we have the following result: 

Proposition 1: If there are matrices P > 0, Ai = diag(An,-- - , Ai„) > 0, A2 = diag(A2i,--- , X 2n ) > 0, 
Q £ R nxn , U as defined above, and a positive real constant p such that the following matrix inequalities hold 

(1,1) PB l + K 1 A 1 -PB 2 +K 2 k 2 

BfP + XiAi -2Ai 0. 

-BlP + K 2 A 2 -2A 2 

(UG ® PD + U ®Q) T + {UG <g> PD + U ® Q) < 0, 
U®P<pI, 

where (1,1) = PA + A T P-Q-Q T + £/, ^ = diag(fc n , • • • ,k ln ),K 2 = diag(fc 2 i, • ■ • , k 2n )- Then the network 
(7) is stochastic synchronization with disturbance attenuation 7. 

Proposition 1 can be proved by replacing F in S2 of (3) by the dynamics of (5), and using the sector conditions 
(6). The details are given in Appendix D. If we choose U beforehand, the matrix inequalities in (8) are all LMIs, 
which are very easy to be verified numerically [15] . For some special G and D, we can further simplify the 
verification process [2TJ [7] . 

An Example 

To demonstrate the effectiveness of our theoretical results, we consider a population of N coupled biological 
clocks, and the individual genetic oscillator is the repressilator [2S]. The repressilator is a network of three 
genes, the products of which inhibit the transcription of each other in a cyclic way (10). Specifically, the gene 
lad expresses protein LacI, which inhibits transcription of the gene tetR. The protein product TetR, inhibits 
transcription of the gene cl, the protein product CI of which in turn inhibits expression of lad, thus forming a 
negative feedback cycle. 

The quorum-sensing system is used for the coupling purpose, which was described in [5]. The system achieves 
cell-to-cell communication through a mechanism that makes use of two proteins, the first one of which (Luxl), 



(8) 
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under the control of the repressilator protein LacI, synthesizes a small molecule known as an autoinducer (AI), 
which can diffuse freely through the cell membrane. When a second protein (LuxR) binds to this molecule, 
the resulting complex activates the transcription lad, as shown in Fig. 1. The noise perturbations in the 
model can arise both intracellular ly, due to the intrinsically noisy property of the gene regulation process, and 
extracellular ly, due to environment fluctuations. 




Figure 1: Schematic representation of the coupled repressilator network. In the left big circle, detailed regulation 
and coupling mechanism are presented. The repressilator module is located at the left of the vertical dotted 
line, and the coupling module appears at the right. 



To model the system, we use a,-, hi, Ci and Ai,Bi, Ci to represent the dimensionless concentrations of the genes 
tetR, ci, lad and their product proteins TetR, CI, LacI, respectively. As in [5] , assuming equal lifetimes of the 
TetR and Luxl proteins, their dynamics are identical, and hence we can use the same variable to describe both 
protein concentrations. The concentration of AI inside the ith cell is denoted by Si. Consequently, the mRNA 
and protein dynamics in the ith cell can be described by [5]: 



da , 
dt 



■didi 

dcj j , a 

dt ~ "3°» fi+B™ 

-diAi + fiiCli, 

-d 5 Bi + p 2 h, 

dSj 
dt 



kSj 



dAj 
dt 

dBj 
dt 

dd 



(9) 



-k s0 Si + k sl Ai + r)(S e - Si), 



where m is the Hill coefficient, S e denotes the extracellular AI concentration, and the meaning of the other 
parameters are standard in genetic network models. We assume that the release of the AI is fast with respect to 
the timescale of the oscillators and becomes approximately homogeneous to establish an average AI level outside 
the cells. In the quasi-steady-state approximation, the extracellular AI concentration can be approximated by 

m 

Qo N 

S e = QpS = — 22 Sj, 
3=1 
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where < Qo < 1 is a constant. Thus the dynamics of Si can be rewritten as 

N 

^ = -{k s0 + (1 - Q )r)}Si + k sl Ai + J2(S j - Si) (10) 

Clearly, the individual model in (9) is of the form (5), in which / = [0, 0, 0, 0, 0, 0, Si/(fi s + Si)) T , Bx is a 7 x 7 
matrix with all zero entries except for B x (3, 7) = k, g = [0, 0, 0, A? /{ix + Af), B^/ifi + B^ 1 ), C™/(^ + C™), 0] T , 
i?2 is a 7 x 7 matrix with all zero entries except for .62(1,6) = #2(2, 4) = 62(3, 5) = a//x, and all the other 
terms are in linear form. Obviously, the coupling term can also be written into the form defined previously. 

The purpose of this example is to demonstrate the effectiveness and correctness of the theoretical result, 
instead of mimicking the real biological clock system. We consider a small size of network with N = 10 
coupled oscillators. The parameters are set as m = 4, a = 1.8, dx — cfe = e?3 = 0.4, // = 1.3, k = 5, fi s = 5, 
d 4 = d 5 = d 6 = 0.5, fa = fa = (3 3 = 0.2,fc s0 = 0.016, k s x = 0.018, Q = 0.8 and r] = 0.4. Since Proposition 
1 holds no matter what Vi(t) is and no matter where it is introduced, and the verification of Proposition 1 
is independent of noise intensity Vi, for simplicity, we set Vi — 0.015 as a scalar for all i, and the noise term 
Viiii(t) is added to the first equation in (9), where rii(t) is a scalar Gaussian white noise process. According 
to Proposition 1 (by letting U — —G, and using MATLAB LMI Toolbox), we know that the above all-to-all 
coupled network can achieve stochastic synchronization with disturbance attenuation 7 = 6. Although 7 is a 
large value, it is easy to show from (2) that the time average of EQIiHj \xi(t) — Xj(t)\ 2 ) is still rather small 
because vf is very small. We omit the computational details here. In Fig. 1 (a)& (b), when starting from the 
same initial values, we plot the time evolution of the mRNA concentrations of tetR (a^) of all the oscillators, 
which behaviors are similar to the experimental results (see, e.g. [1]). Fig. 1 (c) shows the synchronization 
error of a\ — ax for i = 2, • • • , 10. 

In Definition 1, it requires that all the genetic oscillators have the same initial conditions, so that V(x(0)) = 0. 
If the genetic oscillators have different initial conditions, V(x(0)) ^ 0, and thus (12) in the Appendix C is 
replaced by 

J(t) < E{fi[aLV(x{s)) + E^-OfcOMs) " ^(s)) T (xi(s) ~ x 3 (s)) - 7 E< 
+E(V(x(0)). 

Since in genetic networks, the variables usually represent the concentrations of mRNAs, proteins and chemical 
complexes, which are of (not so large) limited values, and so is V(x(0)). For a long time scale, the last term 
of the above inequality is usually much smaller than the absolute value of the first term in the right-hand side, 
and thus the last term can be ignored roughly. In Fig. 3, we show the same computations as those in Fig. 2 
except that the oscillators are with different initial values (randomly in the range (0, 1)). After a period of 
evolution, the network behaviors are similar to those in Fig. 2, which verifies our above argument. In other 
words, rigorously, according to Definition 1, we need that all the oscillators have the same initial conditions, 
but practically, for oscillators with different initial conditions, we can obtain almost the same results. 

For the purpose of comparison, in Fig. 4 we show the simulation results of a network without noise pertur- 
bations. As we can conclude from Figs. 2-4, the networks with noise perturbation, though can't achieve perfect 
synchronization, can indeed achieve synchronization with small error fluctuation, and the network behaviors 
are similar to those of networks without noise perturbations. 
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Figure 2: Simulation results of the coupled repressilators with the same initial values, (a) The evolution 
dynamics of the mRNA concentrations of tetR (cij) of all the genetic oscillators, (b) Zooming in the range 
t £ [600, 700] of (a), (c) The evolution of the synchronization error of at — a\ for i = 2, • • • , 10. 

Synthesis 

In addition to providing a sufficient condition for the stochastic synchronization, Proposition 1 can also be used 
for designing genetic oscillator networks, which is a byproduct of the main results. From the synthetic biology 
viewpoint, to minimize the influence of the noises (on the synchronization), we can design genetic oscillator 
networks according to the following rule: 

min 7, such that the LMIs (8) hold, (11) 

which is obviously from the above theoretical result. This is similar to the synthesis problem in control 
theory. 

Conclusion and Outlook 

In this paper, we presented a general theoretical method for analyzing the stochastic synchronization of coupled 
genetic oscillators based on systems biology approach. By taking the specific structure of genetic systems into 
account, a sufficient condition for the stochastic synchronization was derived based on LMI formalism, which can 
be easily verified numerically. Although the method and results are presented for genetic oscillator networks, it 
is also applicable to other dynamical systems. In coupled genetic oscillator networks, since there is a maximal 
activity of fully active promoters, it is more realistic to consider a Michaelis-Menten form of the coupling terms. 
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Figure 3: Coupled genetic oscillators with different initial conditions: (a) The evolution dynamics of the mRNA 
concentrations of tetR (oj) of all the genetic oscillators, (b) Zooming in the range t € [600,700] of (a), (c) The 
evolution of the synchronization error of <n — a\ for i = 2, • • • , 10. 

As argued in [TJ, our theoretical method is also applicable to this case. To make the theoretical method more 
understandable and to avoid unnecessarily complicated notation, we discussed only on some simplified forms of 
the genetic oscillators, but more general cases regarding this topic can be studied in a similar way. For example: 
(I) The genetic oscillator model (5) can be generalized to more general case such that fi,gi, the component 
of f(y(t)),g(y(t)), are functions of y(t), not only of yi(t), and / and g can also be of non-Hill form, provided 
that < ^jr^jpp < hi,0 < g ^f^p < k 2i ,Va,b e R n (a ^ 6), i= 1, • • • , n, where c u ,c 2i € R n are real 
vectors. (II) Biologically, the genetic oscillators are usually nonidentical. We can consider genetic networks 
with both parametric mismatches and stochastic perturbations in similar ways as those presented in this paper 
and [7J. (Ill) There are significant time delays in the gene regulation, due to the slow processes of transcription, 
translation and translocation. Our result can be easily extended to the case that there are delays both in the 
coupling and the individual genetic oscillators. 

As we know, noises can play both beneficial and harmful roles (for synchronization) in biological systems. 
For the latter case, the noise is a kind of perturbation, and it is interesting to study the robustness of the 
synchronization with respect to noise. In this paper, we addressed this question. For the former case, in 
[131 1 1 4] , the authors studied the mechanisms of noise- induced synchronization. 
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Figure 4: Coupled genetic oscillators without noise perturbation: (a) The evolution dynamics of the mRNA 
concentrations of tetR (oj) of all the genetic oscillators, (b) Zooming in the range t € [600,700] of (a), (c) The 
evolution of the synchronization error of a, — a\ for i = 2, • • • , 10. 

Method 

To simulate the stocahstic differential equaiton x(t) = f(x) + g(x)£(t), the well-known Euler-Maruyama scheme 
is most frequently used, which is also used in this paper. In this scheme, the numerical trajectory is generated 
by + hf(x n ) + Vhg(x n )rj n , where h is the time step and rj n is a discrete time Gaussian white noise 

with < rj n >= and < r) n rj m >= S nm . For more details, see e.g. [18] . 
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Appendices 



A. Notations: 

Throughout this paper, A T denotes the transpose of a square matrix A. The notation M > (<) is used to 
define a real symmetric positive definite (negative definite) matrix. R m denotes the m-dimensional Euclidean 
space; and R nxm denotes the set of all n x m real matrices. In this paper, if not explicitly stated, matrices 
are assumed to have compatible dimensions. E(-) denotes the expectation operator; L2[0,oo) is the space of 
square- integrable vector functions over [0, oo); | • | stands for the Euclidean vector norm, and || • H2 stands for 
the usual L2[0,oo) norm. The Kronecker product A £g) B of an n x m matrix A and a p x q matrix B is the 
np x mq matrix defined as 



A® B 



A n B ■ ■ ■ A lm B 
A n \B ■ ■ ■ A nm B 



For a general stochastic systems 

dx = f(t, x(t))dt + g(x(t))dw(t), 
the diffusion operator L acting on Y(t, x(t)) is defined by 

LY(t, x(t)) = Y t (t, x(t)) + Y x (t, x(t))f(t, x(t)) + -tr&ce[g(x(t))g T (x(t))Y xx (t, x(t))}. 

B. Mean-square synchronization 

For network (1), a natural attempt is to study the mean-square asymptotic synchronization. Analogue to the 
definition of mean-square stability [19j . we can define the mean-square synchronization as follows: 

Definition Al: The network (1) is said to be mean-square synchronous if for every e > 0, there is a 6(e) > 0, 
such that E|xi(£) — Xj(t)\^ >0 < e for \xi(0) — x 3 (0)\ < 5(e), Vi, j. If in addition, lim^ooEla;^) — Xj(t)\ 2 — for 
all initial conditions, the network is said to be mean square asymptotically synchronous. 

In analyzing the synchronization of the network (1), we use the Lyapunov function V(x(t)) = x T (t)(U ®P)x(t) 
[2"T] . where <g> is the Kronecker product, and x(t) = [x[(i),--- 7 x^(t)] T e R Nnyl . According to [2T], this 
Lyapunov function is equivalent to V(x(t)) = Yli<j{~Uij)(xi(t) — Xj(t)) T P(xi(t) — Xj(t)). 

By Ito's formula [19], we obtain the following stochastic differential along (1) 

dV(x(t)) = LV(x(t))dt + 2x T (t)(U <g> P)v(t)dw{t) 
where v(t) = diag(«i, • • • , vn) € R NnxN ; L is the diffusion operator, and 

LV(x(t)) = 2J2 i<j (-U ij )(x i - x^PiFixi) - F( Xj ) - T(x t - Xj )] 
+2x T (t)(U O P)(G ® D + I <g> T)x(t) 
+tra,ce(v(t)v T (t)(U <X> P)) 

We discuss two special cases of the stochastic terms: 

1. The genetic oscillators are perturbed by the same noise, which can occur in the situation that genetic 
oscillators communicate via a common environment. In this case, Vidwi are the same for all i. We let v = 
[vf, • • • , vJ[] T and dw — dwi. Since U is a matrix with zero row sums and u$ is the same for all i, it is easy 
to show that the last term of LV is zero. Thus if the following conditions hold, we will have ~E[dV(x(t))} = 
E[LV(x(t))dt] < 0. 

(Vi - y 2 ) T P[F( yi ) - F{y 2 ) - T( Vl - y 2 )} < 0, 
Vyi,2/ 2 ei?" (2/17^2/2) 
(U <g> P)(G ® D + I ® T) + (G ® D + I ® T) T {U ®P)<Q. 
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Hence, if there are matrices P > 0, T G and J7, such that the above conditions hold, the network (1) will 

achieve mean-square asymptotically synchronization. In this case, roughly speaking, the noise will not affect 
the synchronous state (since they are common for all oscillators), but it will affect the individual oscillator 
dynamics. 

2. The noise intensity matrix Wj is a function of ■ GijDxj, which means that if there is no coupling from 
oscillator j to i, then j does not have contribution to the perturbation of oscillator i. We further assume that 
Vi can be estimated by 



vjvi < 



GijD Xj (t) 



H, 



E 

3 



E 

3 



GijDxjit) 



, ^ > 0. 



Defining v = diag(i>i, • • • , vn) G R 
ing U ® P < pi, we have 

tmce(v (t)v T (t)(U <g> P)) 



< Xr. 



,u>N(t)} T , and H = diag(i?i,- • • ,-ff/v), and assum- 



,{U® P)tr&ce(v(t)v T (t)) 



< P Z>f(*M(*) 



< px T {t)(G ® D) T H{G ® D)a(i). 
So, the conditions for the mean-square asymptotically synchronization of the network (1) in this case are 

(Vi - 2/ 2 ) T P[F(j/ 1 ) - F(y 2 ) - T( yi - j/ 2 )] < 0, 
Vyi,y2&R n {yi^yz) 
(U ® P)(G ® D + I ®T) + (G ® D + I (g>T) T {U ® P) + p(G ® D) T H(G ®D)<0, 

U®P<pI. 

If we consider genetic oscillators of the form of (5), the conditions for the mean-square asymptotically syn- 
chronization can be analyzed by the same method as that in the following Appendix D. 

From Definition Al, we know that the definition of the mean-square asymptotically synchronization is rather 
restrictive, which requires that \vait—>oo^\xi{t) — Xj(t)\ 2 = 0,\fi, j. If it is neither of the above two cases, that 
is neither Vidwi are the same for all i, nor Vi (Vi) reduce to zero when x± = •• • = x n , the network is hardly to 
achieve mean-square asymptotically synchronization. Experimental results also show that usually the genetic 
oscillators can not achieve mean-square synchronization (see for example pQ). So, we argue that the study of 
mean-square synchronization is unrealistic (and therefore meaningless) in genetic networks. 

In Ref. [28] . the authors studied the mean-square asymptotic synchronization of two master-slave coupled 
Chua's circuits. They assume that the noise intensity depends on the difference of the states of the two systems, 
which is also somewhat unrealistic. 

C. Analysis of the general synchronization condition 

To obtain the general synchronization condition (3) of the network (1), we also use the Lyapunov func- 
tion V{x(t)) = x T (t)(U (g> P)x(t). By Ito's formula [19], we obtain the stochastic differential dV{x{t)) = 
LV(x(t))dt + 2x T (t)(U P)v(t)dw(t). According to Definition 1, we assume that the oscillators have the same 
initial conditions, thus we can derive E(V{x[t)) = E (j* LV{x(s))ds^j . For 7 > 0, we define 

J(t) = E{ / [J2(-Uij)(xi(s) - x ] (s)) T (x l (s) - Xj (s)) - 7 E^ T ( s )^( s )] rfs i ( 12 ) 
Then, it is easy to show that for a > 0, 



J(t) < E{f*[aLV(x(s)) + Ei<i(-^i)(^W " 
■(xi(s) -Xj(s)) -7£i v? (s)vi(s)]ds} 
= E{/ * SiOOds}. 



(13) 
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Assuming U <g> P < pi, and letting a = j/p, we have 



Si= ^{2Ei<,(-^)(^W-^(*)) T 

■PlfMt)) - F^t)) - T(a^) - Xj (t))\ 
+2x T (t) {U®P)(G®D + I® T)x{t) 
+trace(v(t)v T (t){U (g> P)) 
+*E i<j (-U ij )(x i (t)-x j (t)) T (x i (t)-x j (t)) 

< z p {E i<j (-U ij )[2(x i (t)-x j {t)) T 
■PiFixJ-FixA-Tixi-Xj)) 
+^{x i {t)-x j {t)) T {x i {t)-x i {t))] 
+2x T (t)(U ® P)(G ® D + I ® I>(t)}. 

If E(5i) < 0, we will have J(t) < 0, and thus, (2) follows immediately from (3). 
D. Proof of Proposition 1 

Proposition 1 can be proved by replacing F in S 2 of (3) by the dynamics of (5), and using the sector conditions 
(6). We have 

S 2 = 2{ yi {t) - y 2 {t)) T P[{A T){ yi {t) y 2 {t)) 

+B 1 (f(y 1 (t)) f(y 2 (t))) B 2 (g( yi (t)) g(y 2 (t)))} 

+^yi(t)-y 2 (t)) T ( yi (t)-y 2 (t)) 
< 2( yi (t) - y 2 (t)) T P(A T)( yi (t) y 2 (t)) 
+£(yi(t)-y 2 (t)) T ( yi (t)-y 2 (t)) 

+2( yi (t) lto(t)) r PBi(/(i/i(*)) - f(y 2 (t))) 
-2( yi (t) - y 2 (t)) T PB 2 (g( yi (t)) - g(y 2 (t))) 

-lEUMMvuit)) - Myn(t))) 
■[(fi(yu(t)) - fiivnit))) - hiyuit) - y 2l (t))} 
-2Y,U^ii9i{yu{t))-gi{y2i{t))) 

■[(9l(yil(t)) ~ 9i(y2i(t))) - k 2 {y u {t) - y 2l {t))]. 
By letting Q = PT and denoting the first matrix in (8) by M\, we have S 2 < £(t)Mi£(t) < for all yi,y 2 
except for yi = y 2 , where = [( Vl (t) - y 2 (t)) T , (f(yi(t)) - /(y 2 (t))) T , (ff(i/i(t)) - g(y2(t))) T } T G R 3nxl . So, 
the first condition in (3) is satisfied. Substituting Q = PT, the second inequality in (3) is equivalent to the 
second inequality in (8). Thus, Proposition 1 is proved. 
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